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Abstract 

The effect of thermally generated bulk stochastic forces on the statistical 
growth dynamics of forwards bifurcating propagating macroscopic patterns is 
compared with the influence of fluctuations at the boundary of a semiinfinite 
system, < x. To that end the linear complex Ginzburg-Landau amplitude 
equation with additive stochastic forcing is solved by a spatial Laplace trans- 
formation in the presence of arbitrary boundary conditions for the fluctuations 
of the pattern amplitude at a; = 0. A situation where the latter are advected 
with an imposed through-flow from an outside upstream part towards the 
inlet boundary at a; = is investigated in more detail. The spatiotemporal 
growth behavior in the convectively unstable regime is compared with recent 
work by Swift, Babcock, and Hohenberg where a special boundary condition 
is imposed. 
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I. INTRODUCTION 



Nonequilibrium extended systems that undergo a pattern forming instability with 
nonzero group velocity — whether the latter arises in open-flow configurations via an 
imposed through-flow @-|l5| or internally from a bifurcation to travelling waves [[l6l -|T8| — 



show a peculiar sensitivity to perturbations in the so called convectively unstable parameter 
regime . Therein initially spatially localized perturbations superimposed on some ba- 

sic unstable, macroscopically unstructured nonequilibrium state of the system are advected 
downstream by the group velocity. Simultaneously they are selectively amplified [^,^ ac- 
cording to the deterministic amplification mechanism that occurs near the pattern forming 
bifurcation Thus, advection spatially deconvolutes the pattern growth from its cause - 
the perturbation source. But the latter leaves its signature via the bifurcation-induced, de- 
terministic amplification process on the spatiotemporal properties of the downstream grow- 
ing macroscopic nonequilibrium structures. The statistical dynamics of fluctuations of these 
structures is easily detectable p|-3,p!0|- [12| , p!7| , pr8| thus providing information about the small 
perturbations that have triggered their growth. 

The feasibility to distinguish the influence of different perturbation sources - e.g. those 
that are external to the system under study versus internal thermal noise - on the macro- 
scopic pattern growth in a very direct way has attracted much research activity recently 



,pI) Hl^ , M ,[rq] • Here one has to keep in mind that the deterministic bifurcation dynamics 
of the macroscopic structure forming instability determines which of the modes that a per- 
turbation source emits is amplified and which one decays in downstream direction. And as 
long as the mode amplitudes are still small each of the perturbation modes evolves separately 
and independently: the bifurcation dynamics is initially linear with saturating nonlinearities 
being irrelevant thus producing exponential growth with time and downstream distance. 

Hence, in a system with group velocity in positive x direction one expects that fluctu- 
ations generated near the inlet boundary at x = and/or perturbations that might enter 
the system via the inlet are most important. They have the longest advection/amplification 
time to grow in the flow before they reach a particular observation point at a downstream 
location x > 0. Therefore a careful analysis of the influence of the inlet boundary condition 
relative to the effect of bulk (fluctuating) forces is necessary in particular in view of the fact 
that both sources cause exponential growth. Such an analysis should be able to shed some 
light on the problem whether in bounded experimental systems it is internal bulk gener- 
ated thermal noise that triggers and sustains the growth of macroscopic vortex structures in 
the convectively unstable regime or other kind of environment-produced perturbations that 
might, e. g., be advected from the inlet boundary into the bulk. 

On the theoretical side the effect of noise on the pattern forming process has been inves- 
tigated mainly within the framework of the complex Ginzburg-Landau amplitude equation 
(CGLE). It describes the slow spatio-temporal behavior of the complex amplitude A{x,t) 
of the critical mode of the pattern that can grow first when crossing the bifurcation thresh- 
old. Bulk thermal noise is incorporated in this approximation scheme by fluctuating volume 
forces f{x,t) that act additively on A{x,t) and that are derived p^-p6| for local thermal 
equilibrium conditions of driven systems from the linear Landau-Lifshitz equations of fluc- 
tuating hydrodynamics around global thermal equilibrium |2^. But the problem of the inlet 
boundary condition that arises with the CGLE in a finite domain has been addressed only in 
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some idealized way: Deissler p8| imposes A = ai the inlet. Liicke and Recktenwald [p9[ 
vestigate the pattern growth resulting from an inlet boundary condition that is produced by 
thermal equilibrium transverse momentum fluctuations which are advected into the system 
at X = 0. Babcock, Ahlers, and Cannell have solved the stochastic CGLE numerically in 
a system that was augmented in upstream direction, i.e., to negative x by a small subcritical 
part and compared with imposing A = at x = 0. Swift, Babcock, and Hohenberg (SBH) 
23| consider solutions of the stochastic CGLE where the amplitude fluctuations A{x = 0,t) 



at the boundary are strongly correlated with the bulk thermal forcing /. 

In this paper we first of all solve the stochastic CGLE subject to arbitrary inlet boundary 
conditions with a spatial Laplace transformation method. We investigate the influence 
of thermally generated bulk additive forcing on the statistical dynamics of downstream 
growing amplitude fluctuations in comparison with the effect of statistically stationary inlet 
fluctuations that are independent of the bulk forcing. Furthermore, we use a simple model to 
estimate fluctuations that are swept with an imposed through-flow into some experimental 
systems via the inlet. This allows to assess more quantitatively the importance of such inlet 
fluctuations relative to bulk thermal forcing. Experimental setups to tune and to control 
the former are shortly discussed as well. Finally we compare in detail with the results of 
SBH. 

Our paper is organized as follows: In Sec. the linear stochastic CGLE with thermally 
generated bulk additive forcing is briefly reviewed. The solution for the fluctuating ampli- 
tudes in frequency space, A{x,uj), and their correlation spectra are given in Sec. |T|. The 
next section presents the model for inlet fluctuations and compares their downstream evo- 
lution with the growth of bulk generated fluctuations. The last section contains a summary. 
In Appendix ^ we present the " space- retarded, causal" Green function of our solution in the 
convectively unstable regime. In Appendix ^ we compare our results for Green function, 
amplitude fluctuations, and correlation spectra with those of SBH. 



II. SYSTEMS 

The systems that we have in mind show a continuous, i.e., hysteresis- free transition, say, 
from a spatially homogeneous state to a spatially structured one. Thus, a spatially periodic 
solution branches in a forwards bifurcation out of the basic solution at a critical control pa- 
rameter. The most intensively studied examples for such systems are the Rayleigh-Benard 
(RB) problem of straight roll convection vortices in a fluid layer heated from below |T[] or 
the rotating Taylor-Couette (TC) problem of axisymmetric toroidal Taylor vortices in the 
annulus between concentric cylinders of which the inner one is rotating [|T],0. But our in- 
vestigation being based on the CGLE for such ID patterns also includes other problems for 
which this equation captures the essential parts of the pattern formation process. Neverthe- 
less we shall use in the remainder of this paper the terminology of the RB or TC system. 
Thus, we consider in particular open RB (TC) setups with an externally imposed horizontal 
(axial) mean flow through a not too wide rectangular convective channel (annulus between 
concentric cylinders). Then the deterministic field equations show an oscillatory instabil- 
ity. The structureless, stationary basic flow without vortices is stable for Rayleigh numbers 
Ra (Taylor numbers T) below a critical threshold RadRe) [Tc{Re)] that depends on the 
through-flow Reynolds number Re. At the critical threshold a stable, spatially extended. 
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temporally oscillatory state of downstream propagating vortices branches off the basic state. 
The latter is unstable and the former is stable for supercritical driving. 
We use the relative control parameter 

Ra T 

to measure the distance from the onset of propagating vortex flow for Re ^ and of 
stationary vortex flow for Re = 0, respectively. In this notation fic = 0, i.e., 

-ROf. (-Re) , , Tf. (-Re) 

is the critical reduced threshold for onset of propagating vortex flow. The shear forces of the 
through-flow slightly stabilize the homogeneous basic state, so ec (-Re) slightly increases with 
-Re PJ^J^-|8|,|3T1| . Note that in the RB system we consider so-called transverse vortices whose 
roll axes are enforced to lie perpendicular to the through-flow direction by the geometrical 
design of the rectangular convection channel. So for not too large -Re the critical wave vector 
of the vortex patterns is oriented in both cases along the x-direction of the through-flow. 



A. Amplitude equation approximation 

For small supercritical control parameters, < ;U ^ 1, the macroscopic vortex dynamics 
is governed by a narrow band of near critical modes. This allows to represent the hydrody- 
namic fields within an amplitude equation approximation by a complex amplitude A{x, t) 
multiplying the critical mode 

w {x, rr, t)=A {x, t) e^C^-^-^-*)^!; (r^) + c.c. (2.3) 

Here w is for example a component of the vortex velocity field. The critical wave number kc, 
frequency Uc P,|6|-p|,pll-^ , and eigenfunction w{r±) [p2| , p4| depending on the coordinates 
perpendicular to the critical wavevector of the downstream propagating vortex structures 
are known as functions of -Re from a linear stability analysis of the basic flow. 

In the absence of fluctuations the spatiotemporal behavior of the complex vortex ampli- 
tude A{x,t) is determined by the ID CGLE 

'To{dt + vdx)A(x, t) = (1 + ico)/i + (1 + icij^^dl. A(x, t) + nonlinear terms. (2.4) 

The coefficients tq and are even in -Re while the group velocity v and the imaginary parts 
Co, Ci are odd in -Re P,pT[]. 

B. Absolute and convective instability 



For small /i and -Re the control parameter plane is divided into three domains char- 
acterized by different growth behavior of linear perturbations of the basic flow state. In 
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the presence of through-flow one has to distinguish [|T9[ between the spatiotemporal growth 



behavior of spatially localized perturbations and of spatially extended ones. Below the 
threshold /i = for onset of propagating vortex patterns any perturbation, spatially local- 
ized as well as extended, decays. This is the parameter regime of absolute stability of the 
basic state. 

For yU > extended perturbations can grow. However, a spatially localized perturbation, 
i.e., a packet of plane wave perturbations is advected in the so-called convectively unstable 
parameter regime faster downstream than it grows — while growing in the comoving frame 
it moves out of the system [^^-|2T],^ . Thus, the downstream as well as the upstream facing 
intensity front of the vortex packet move into the same direction, namely, downstream. A 
spatially localized perturbation is blown out of any system of finite length and the basic 
state is reestablished. It therefore requires a persistent source of perturbations like, e. g., 
noise operating in the bulk or fluctuations entering the system via the inlet to sustain a 
vortex pattern in the convectively unstable regime [^. The spatiotemporal behavior of 
such " noise sustained structures" P,|l2| being determined by the noise source thus allows to 
infer properties of the latter from the former. 

In the absolutely unstable regime a localized perturbation grows not only in the down- 
stream direction but it grows and spatially expands also in the upstream direction until the 
upstream propagating front encounters the inlet in a finite system The final nonlinear 



stucture i,|,|,|T0|JT3-|T|,||] resulting in such a situation is self sustained and stable. For 
sufficiently large amplitudes this structure is dominated by the deterministic contributions 
in the relevant balance equations and thus it is insensitive to, say, thermal noise. 

Within the framework of the amplitude equation the boundary between absolute and 



convective instability is given by ||2T[| 



^^^'^ " 4(l + cf)eo'- ^^'^^ 

Thus, the convectively unstable regime investigated here is characterized by < /i < fi^onv 
or, equivalently, by the reduced group velocity 

V = v = 2J^™^ (2.6) 




being larger than 2. It should be mentioned that the amplitude equation approximation 
of fJ^conv describes the boundary between absolute and convective instability resulting 



from the full hydrodynamic field equations very well in the RB and TC systems for 

the small Reynolds numbers considered here. 



C. Stochastic amplitude equation 

The starting point of most theoretical approaches for describing the effect of thermal 
noise on macroscopic hydrodynamic pattern growth is the Langevin concept of Landau and 
Lifshitz ||2^. Therein the linear hydrodynamic field equations are heuristically supplemented 
by Gaussian distributed statistically stationary stochastic forces with zero mean. Their two- 
point correlation functions refiect the isotropy and translational invariance of an unbounded 
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system and their frequency and wave vector spectra are white. Projection p2|-p6[] of these 
equations onto the critical modes leads then to the linear stochastic CGLE 



Toidt + vd,) - (1 + zco)/i - (1 + tc,)eod: A{x, t) = fix, t) 



(2.7) 



with Gaussian distributed complex forces f{x,t) that are 5-correlated in space and time 



<f{x,t) >=< f{x,t)f{x',t') >=0, 
< fix, t) [fix', t')T >= F5ix - x')5it - t'). 



(2.8a) 
(2.8b) 



Here the star implies complex conjugation. The noise strength F has been evaluated by 
several authors p^-p6| for the RB and TC systems without through-flow and with different 
boundary conditions. Since we do not need the specific values of F we shall not dwell on 
questions whether and how the forces are changed by finite-size effects, system geometry, and 
the through-flow. We rather take eqs. (|2.7| , |2.8|) as the starting point for our investigations. 

However, for later use we mention that our F is related via F = 2C,oTqF^^^ to the 
noise strength F^^^ given in eqs. (2.10 - 2.16) of SBH. We should also like to mention 
that the value of F characterizing a particular system depends on the arbitrary choice of 
the normalization condition for the linear critical mode w (r^) that enters into Eq. ( p.3|) 
and that establishes the relation between the physically observable order parameter field 
w ix, r±; t) and the complex amplitude A of the critical mode. 



III. LINEAR STOCHASTIC CGLE IN SEMIINFINITE SYSTEMS 

In this section we discuss the statistical dynamics of pattern growth in the convectively 
unstable parameter regime as described by the stochastically forced linear CGLE (2/7). We 
consider a semiinfinite system with the through-flow entering the system at x = 0. This 
flow transports fluctuations into the system that we specify by imposing an inlet boundary 
condition A (x = 0, t) on the amplitude A. We consider situations in which A (x = 0, t) is the 



result of a statistically stationary process; see Sec. |IV A| for a model. Thus, with the bulk 
forces fix, t) being also statistically stationary we restrict our investigations to the long-time 
statistically stationary fluctuations of Aix,t) — initial transients occurring in experiments 
after, e.g. changing control parameters have died out. Under these circumstances we found 
it convenient to perform part of our analysis in frequency space. After a temporal Fourier 
transformation 



/oo 
dte'^^'Aix, t) 
-oo 



(3.1) 



the stochastic CGLE becomes an ordinary, second order, inhomogeneous differential equation 



Toi-iuj + vd^) - (1 + ico)^i - (1 + ici)^od^ A(x, uj) = fix, u) 



(3.2a) 



that is to be solved subject to the inlet boundary condition y4(x = 0, uj) 
derivative of A at the inlet or, equivalently, on the behavior of A for x 
a physical requirement discussed in Sec. [HI A| below. 



A condition on the 
+00 follows from 
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For later use we rwrite the CGLE in the form 

- (1 + tc,)eo{d. - tK,){d^ - tK2)A{x,uj) = f{x,uj) (3.2b) 

that displays the spatial eigenvalues K12 of the CGLE explicitly — in the solution of Eq. 
( p.2|) there appear contributions of the form e^^^^ and e*^^^ where Ki^2 are the complex 
a;-dependent characteristic exponents, i.e., spatial eigenvalues of Eq. 



^0 VI + iciKi = ±i^J (1 - ici)n%^^ - (1 + «co)/U - iutq - - ici)/i^„„^. (3.3) 

The physically relevant characteristic exponent that controls the spatial growth of vortices 
in downstream direction is Ki whereas the second eigenvalue of the CGLE is such that 
ImK2 < for all fi and all u. 

ImKi{^;uj) is for < positive so that the spatial growth rate —ImKi is negative 
implying spatial decay of each frequency mode uj. For /i > selective amplification of 
frequency modes out of a finite band < u < uj^ occurs 

ImKiioj; < fi < jJi^conv) < for c<j_ < c<j < c<j_|_ (3.4a) 
ImKiiuj; < fi < ^l^^^) > for u < and u > uj+. (3.4b) 

The band limits for amplification are 



ToUJ± = ±2^ (1 + cf )/i/i^„„^ + (ci - co)/i. (3.5) 

ImKi has a quadratic minimum at the frequency 

= -(co + Ci)yu/ro (3.6) 

which is, due to the smallness of cq and ci, in general located close to u = 0. Thus, the 
frequency mode with largest spatial growth is ujm- 

For plots of —ImKi versus u see, e.g., [^, Fig. 1] and Fig. 14]. The relations to our 
exponent Ki are iKi = k [M = p M. 



A. Spatial Laplace transformation 

We have solved the CGLE ( ^.2| ) with the method of Green's functions and, equivalently, 
using a spatial Laplace transformation. The latter method most clearly displays the influence 
of the inlet boundary properties on the solution. Thus, we start with the latter method here 
and present the relation to the former in Appendix 

We use the spatial Laplace transformation 

A{K) = -I \ dxe-'^'^Aix) ; K complex. (3.7) 
Jo 

Here the second argument - u; or t - of A is suppressed. Using the transformation properties 



d,A{x) ^ iKA{K) + zA(x)U=o, 

dlA{x) ^ -K^A{K) + i{d, + tK)A{x)\,=^ 



(3.8a) 
(3.8b) 
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the Laplace transformed CGLE (|3.2|) reads 

(1 + tCi)eo{K - K^){K - K,)A{K, u) = 
i{l + tci)eo [dx + i{K -K,- K2)] A{x, uj) U=o +fiK, lu). (3.9) 

Its mathematically formal, physically unrestricted solution 

z[d.^ + z{K-K,-K,)]A{x,uj) U=o f{K,u;) 

^ ' ^" (K-K^KK-K,) + {l + ^c,)a{K-K,){K-K,) ^'-'"^ 

contains poles in the complex i^-plane at Ki and at K2. 

Now in the physical solution the residue of the pole at K2 has to vanish, i.e., 

hni [iK-K2)AiK,uj)] = 0. (3.11) 

Otherwise there would be a contribution ~ g*^23; f^j, _^ giving rise to spatial growth 
and a divergence of A{x +00,0;) even for subcritical driving /i < 0. The physically moti- 
vated requirement that there is no growth ~ e*^^^ for x — » +00 is equivalent to selecting 
the "retarded" or causal Green function in the solution of the CGLE (cf. Appendices ^ and 
|B 1| ). Therein the solution is represented by a superposition of "waves" which propagate 
(in the frame comoving with velocity v) outwards, i.e., away from a localized perturbation 
source ~ 6{x — xo)6(t — to) that creates a perturbation pulse at xq, to. For < the 
"waves" are damped while for < fi < ficonv they grow in downstream direction within the 
frequency band (|3.5|) . On the other hand, the "advanced" Green function describes inwards 
propagating waves ~ e*^^^'^'^^^"'^*! that create the pulse at xq, to- They would cause even 
for /i < a divergence of A at x — > +00. 

Note that the requirement ( |3.11| ) implies that for a given forcing f{x, u) the initial slope 
dxA{x,Lj) of the amplitude at the inlet is fixed by the inlet value A{x = 0,0;) and 
the forcing f{K2,uj) = —i dxe~^^^^ f{x^uj). From eqs. ( |3.10D and ( pj.llp one obtains the 
relation 

t{d,-tK,)A{x,u) 1,^0 = 0. (3.12) 

Using this restriction to eliminate the slope in (|3.10|) the physical solution of the CGLE can 
be written as the sum of two terms 

A{K,uj) = A,n{K,uj) + Ab{K,uj). (3.13a) 

Here 

AUK.u) = -^^^^^ (3.13b) 
is the contribution from the inlet fluctuations A{x = 0,uj) while 
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is caused by bulk forcing. In real space the solution reads 

A{x, uj) = Ain{x, u) + Ab{x, u) 

where 

Ain{x,uj) = A{x = 0,o;)e^^^^ 

and 



Ab{x,u) 



[l + ^c^)eo{K^-K2) Jo 



dx' 



^iKi{x—x') ^iK2(x—x') 



^iKix _ ^iK2X 



f{x\u) 
dx'e-'^''""' f{x\uo). 



Note that the bulk part vanishes at the inlet boundary, Ab{x = OjCj) = 0. 



(3.14a) 
(3.14b) 



(3.14c) 



B. Digression: constant bulk and inlet forcing 

In this subsection we make a digression to investigate the special case 

/(x, uj) = 2n5{uj)fo ; A{x = 0,uj) = 27cS{uj)Ao (3.15) 

of constant (deterministic) forcing f{x,t) = /o in the bulk and constant inlet boundary 
condition A{x = 0,t) = Aq. This investigation is useful and instructive for two reasons: 
(i) It allows to check the general solution ( |3.14| ) for a special case that can easily be solved 
analytically in a direct way without Laplace transformation and that furthermore can easily 
be solved numerically, (ii) It gives insight into the way how initial slope dxA \x=o, inlet value 
Aq, and bulk forcing /o are tied to each other. 

We consider in this subsection the stationary solution A{x) of the CGLE , i.e., A{x, u) = 
2n6{uj)A{x). Then (^A^ gives 

An(x) = Aoe^^?" (3.16a) 

and ( |3.14c| ) becomes 

/o 



Ah(x) 



(l - e^^?^) 



(3.16b) 



with = Ki{uj = 0) and K2 = K2{uj = 0) given by ( |3.3| ). The complete solution reads 
A{x) = Ain{x) + Ah{x) 



Aoe 



iK°x 



(1- 



AK°x 



(3.16c) 



{l + ico)n 

Here we have used (1 + ici)^QKlK2 = —(1 + ico)fi. The physically motivated restriction 
( p.ll| , p.l2| ) on the slope at the inlet becomes in our special case 

/o 



dxA{x) \x=o= iK^Ao - i- 



An + 



fo 



(3.17) 



1 + ico)fi 

since f{K2, uj) = —fo2TT6{uj) / K2 for constant forcing. We have checked the validity of ( p.l6c|) 
also by direct numerical integration of the (nonlinear) CGLE. 
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C. Correlation functions 



Here we investigate the frequency spectrum 



d{t - t')e'"(*^*'^ C{x, \t-t' 

-oo 

of the time-displaced autocorrelation function 

C{x, \t-t'\)=< A{x, t) [A{x, t')]* > 



(3.18) 



(3.19) 



of amplitude fluctuations A{x,t) at the same downstream position x. Since the semiinfi- 
nite system with through-flow is spatially not translational invariant the correlations (|3.19|) 
depend explicitly on the position x. On the other hand, the statistical stationarity of the 
forcing processes in the bulk and at the inlet causes a temporal dependence only via | t — t' \. 
Using this property one has 



C{x, u) 



< A{x,uj) [A{x,uj')Y > 
2n5{uj - to') ' 



(3.20) 



We consider inlet fluctuations A{x = 0,a;) that are determined by forces operating 
outside the bulk fluid volume, for instance by fluctuating forces in the upstream part of 
an experimental setup. Thus, in the convectively unstable regime the contribution Af,{x,uj) 
( ^.14cD due to bulk forcing is uncorrelated with the inlet fluctuations A{x = 0, u) so that 



< A{x = 0,uj) [Ah{x,uj')]* >= =< Ain{x,uj) [Af,ix,uj')]* > 



Therefore 



C{X, U) = Cin{x, U) + Cbix, U) 



(3.21) 



(3.22a) 



is the sum of correlations 



Cin{x,Uj) 



-2ImKix 



D(uj) 



(3.22b) 



of the contribution Ain{x,u!) ( p.l4b|) caused by inlet fluctuations and of correlations 



Cb{x,uj) 

+ ( 4/m 



F r 1 1 

2(1 + c?)e^ \K^-K2\^\ IniK^ ~ ImK^ 
1 1 1 ~ 



-2ImK\x 



+4Re 



Kl - K2 ImKi ImK2 
1 I 



i(Ki-K2)x 



2ImK2 KI-K2) 
of Ab(x, cu) ( |3.14c| ) caused by bulk forcing. Here 

< A{x = 0,uj) [A{x = 0,uj')]* > 



(3.22c) 



D{uj) = Cin{x = 0,u) 



2n5{u - Lu') 



(3.23) 
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is the fluctuation spectrum of the inlet amplitude. The expression ( |3.22c| ) for 



a(...) = li^(£41M£ii2i:ii. (3.24) 

is obtained directly from ( p.l4c| ) by using the force correlations 

< /(x, uj) [fix', uj')]* >= 6{x - x')2tt6{uj - u')F. (3.25) 

Let us briefly discuss the small-x and the large-x behavior of C{x,uj): Since Cb{x = 
0,u) = one has C{x = 0,uj) = Cin{x = 0,u) = D{uj). Furthermore, Cj„ varies linearly 
and Ch quadratically in x close to the inlet. 

The large-x behavior is physically more interesting since it tells how and which frequency 
modes are amplified. Here one has to distinguish the two cases (i) ImKi > and (ii) 
ImKi < 0. The case (i) holds for all u; if < and for uj outside the band {uj-,uj+) of 
growing frequency modes if /i > 0. For these damped modes the autocorrelation approaches 
at large distances from the inlet the limit value 

lim C{x, uj) = , ^ f - 7-^) . (3.26a) 



2(1 + 4)^0 \K1-K2 |2 \ImKi ImK2. 

The fluctuations of these damped modes are influenced and sustained at large distances 
from the inlet only by the bulk forcing. The limiting value ( |3.26aD coincides with the 
autocorrelation 

/oo p 
-oo 2^ I (1 + ici)ilk^ - /^(l + ^Co) - i{uj - vk)To [^ ^^''^^^^ 

of amplitude fluctuations in an infinite system, —00 < x < 00, that is translational invariant 
without inlet boundary conditions. 

In case (ii) of ImKi < 0, i.e., for /i > and uj- < uj < there is exponential growth. 
Then the correlation function is dominated for large x by 

F /II 



C{x ^ 00, u) 



D uj) - 



2(1 + cj)Q I - K2 p \ImKi ImK2 
To derive ( |3.27| ) from ( p.22| ) we have used the relation 



Aim- ^ 



\Ki-K2\^\ Kl - K2 ImKi ImK2j \ Kl - K2 P KlmK^ ' ImK2, 

(3.28) 

Eqs. (|3.22|) for the spectra of amplitude fluctuations and their large-x limiting behavior 
are the central result of this paper. It allows to assess the influence of (externally generated) 
inlet fluctuations and of bulk stochastic forces separately. 



IV. INLET AND BULK GENERATED FLUCTUATIONS 



Here we first introduce a simple model for the frequency spectrum D{uj) ( 3.23 ) of inlet 
fluctuations of the amplitude A{x = 0). It allows to make quantitative estimates of the 
importance of their contributions Cin{x,uj) (p.22b| ) relative to those of the bulk generated 
fluctuations Cb(x, a;) ( p.22c|) . 
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A. Model for inlet fluctuations 



In Ref. |2^ D{uj) was estimated by projecting transverse momentum current fluctuations 
in thermal equilibrium onto the critical mode of the hydrodynamic field equations. These 
fluctuations were assumed to be advected towards the inlet from an upstream part x < 
without deterministic driving, fi = —1, and to enter the main supercritical part of the 
experimental cell , e.g., via a porous plug at x = 0. Thus, the Rayleigh number Ra or 
the Taylor number T in the upstream part was taken to be zero in |^9[. In experiments, 
however, one often has an upstream section — L < x < to the left of the inlet at x = in 
which the fluid is equilibrated at a finite subcritical driving /i < 0. In the remainder of this 
paper we mark quantities characterizing such a subcritical upstream part of the system by 
an underbar. 

Babcock, Ahlers, and Cannell have numerically solved the stochastic CGLE for such 
a setup with L = 0.1, = —0.1, and boundary condition A{x = —L,t) = 0. We consider a 
situation where L is sufficiently large. Then perturbations, e.g., from the pumping machinery 
that are swept at x = —L into the subcritical upstream part of the apparatus have decayed 
by the time the fluid reaches the inlet, x = 0, of the supercritical part. Thus, one ideally 
has at X = only fluctuations in the hydrodynamic fields that are sustained by thermally 
activated fluctuating forces. The most important field fluctuations are those with wave 
numbers close to the critical one and with spatial variation in the directions perpendicular 
to the through-flow similar to the critical mode. We approximately describe the dynamics 
of the amplitude of these modes with the linear stochastic CGLE ( |2.7| ), now however for 
subcritical driving yU < 0. Therein, the spectrum of amplitude fluctuations at x = that 
are generated at subcritical driving by thermal forcing /(x < 0,t) alone — L is large — is 
given by dpI) with /X = /X < 

The total spectral power of these noise generated fluctuations of the amplitude A{x = 0) of 
the critical mode at the inlet 



increases with increasing driving in the subcritical upstream section and diverges in the 
absence of nonlinearities in ( p^.7| ) at the critical driving /i = 0. Thus, the estimate ( [4.1|) is 
only realistic for a situation where the noise generated amplitudes in the upstream part are 
small enough to neglect the nonlinear contributions in the hydrodynamic field equations, 
i.e., for /X not too close to zero. If that is the case then a variation oc 1/y^— /I in the inlet 
contribution should be experimentally observable. 

In Fig. |l|we show the inlet spectrum Diu) (^?T| ) normalized by ToD{t = 0) versus reduced 
frequency utq for several subcritical upstream control parameters as indicated. The inlet 
spectrum D{uj) is dominated by the term 1/ImKi. It controls the /i variation of the peak 
height and of the peak width of D{uj) and it governs the strong enhancement of those modes 
around u = that can grow in the supercritical downstream part of the system. 
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Finally we should like to comment on the objection p5| against using in a system with- 
out translational invariance thermal equilibrium spectra for inlet fluctuations that come 
from the fx = —1 upstream part in |2^. This seems to be well justified in view of the fact 
that downstream forces do not influence upstream locations — downstream generated per- 
turbations cannot propagate upstream if the through-flow is smaller than the convectively 
unstable threshold. Thus for any x in the upstream section that is sufficiently far away from 
L one should observe fluctuations as in an infinite system. For the stochastic CGLE this is 
reflected by the equality of (|3.26a]) and ( p.26b|) . But the corresponding property holds also 
for the full linear field equations with through-flow. 



B. Comparison of inlet and bulk generated fluctuations 

In Fig. D we show the downstream evolution of the spectral contributions Cin{x,uj) 
( ^.22b| ) and Cb{x,u) ( |3.22c| ) from inlet and bulk forcing, respectively, to the total spectrum 



of amplitude fluctuations. The contributions are reduced by F/^q measuring the forcing 
strength. The inlet fluctuations C{x = 0, a;) = Cin{x = 0,u) = D{u) are taken from the 
model of the previous subsection with a subcritical upstream driving of /i = —0.1. The 
Reynolds number is Re=2 and the supercritical downstream driving is = /U^q„^/2 = 0.03. 
For these parameters the band limits (|3.5| ) of amplification ImKi{uj) < are uj-Tq = —0.084 
and uj+tq = 0.086, respectively. While the logarithmic scale of the ordinate tends to obscure 
the boundaries between growing and decaying frequency modes the left border U- of the 
amplification band can be identified in the plots of Cin{x, uj) of Fig. ^ by their crossing point 
— there Cj„-spectra at different x have the same value for ImKi = . 

The bulk contribution Cb which vanishes at the inlet is at a distance of one correlation 
length ^0 still smaller than the inlet contribution. Then, with growing downstream distance 
Cb becomes larger. At x = 20^o the bulk contribution has almost reached its large-x single 
exponential growth asymptotics oc ^-'^^'^^li^)^ ^ There the contribution to Cb from the last 
term oc e^[^i('^)-^2(^)l^ in ([3:22^ ) is already small. 



That can also be seen from Fig. y. There we show for the same supercritical downstream 
driving /i = /i^Q„^/2 = 0.03 that was used in Fig. ^ the downstream growth of the spectral 
peak Cb{x, u = 0) in comparison with the peak height Cin{x, u; = 0) of the inlet contribution 
— for our parameters the spectral maxima are located practically a.t uj = 0. Cm was 
evaluated for several subcritical upstream driving values of /i as indicated. The dashed line 
labelled Csbh is discussed in Appendix [B 3|. Thus, for the values of n in Fig. H the bulk 



contribution rapidly outgrows the inlet contribution. A similar observation was made also 
in the numerical treatment of Babcock, Ahlers, and Cannell [§. 

The slope of the straight lines Cin and of the asymptotics of Cb and Csbh is given by 
—2ImK^ where the superscript indicates that the eigenvalue Ki is taken at = 0. In the 
far downstream region the ratio of the spectral maxima of Cb and Cj„ approaches according 
to ( I3.27D the limiting value 



a(x,^ = 0) _ F/D{u = ^) / 1 1 \ 

™ Cin{x, = 0) 2(1 + cDit I Kl* - Kl |2 [imK^ ImK^ ' ^ ' 



If one uses for D(uj = 0) the result ([4.1|) of our model one finds for the parameters of Fig.|^ 
that the bulk generated downstream fluctuations Cb dominate over the fluctuations that 
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are swept into the system through the inlet whenever the upstream driving < —0.025. 
If the downstream driving /i is smaller than in the example of Fig.|^ then the size of inlet 
fluctuations that would ensure Ctn to be comparable to Cb for large x would have to be 
larger - i.e., would have to be even closer to zero. 

Thus we conclude that typically the large (small)-x behavior is governed by bulk (inlet) 
fluctuations. However, our upstream model clearly shows that in appropriately designed ex- 
perimental upstream sections — with L being long enough to avoid influence of the pumping 
machinery operating at a; = — L — one can tune and enhance the inlet noise at x = almost 
arbitrarily by selecting the subcritical upstream driving n accordingly. Finally we should 
like to stress that our model for inlet fluctuations contains only internal thermal fluctuations 
generated in the subcritical upstream part. It does not contain other perturbations, e.g., 
from the pumping machinery that might in an experimental setup be advected towards the 
inlet. Their contribution to the inlet fluctuation spectrum would have to be added to our 
model spectrum ( [4.1| ). However, since our general formulas in Sec. JTT|have been derived for 
arbitrary inlet fluctuations with spectrum D{uj) they also apply to such a situation. 

V. CONCLUSION 

We have investigated the effect of inlet fluctuations and bulk stochastic forcing on the 
noise sustained pattern growth in semiinfinite systems, < a;, with through-flow in the 
convectively unstable parameter regime. We have used a ID CGLE that is appropriate 
for forwards bifurcating structures in order to quantitatively compare both effects on the 
statistical dynamics of the downstream growth of the complex pattern amplitude A{x,t). 

To that end we have solved the linear stochastic CGLE subject to arbitrary, statistically 
stationary boundary conditions at the inlet x = in the presence of bulk additive forcing 
f{x,t) with a spatial Laplace transform. This method is well suited for the semiinfinite 
geometry: inlet and bulk forcing are easily dealt with separately since they both appear 
explicitly and additively via boundary and inhomogeneous terms, respectively. Furthermore, 
this approach easily allows to identify the contributions to the general mathematical solution 
from a pole at K2{uj) in the complex wave- number plane that have to be discarded on 
physical grounds since they would cause growth even for subcritical control parameters. The 
physically motivated restriction is equivalent to choosing the " space- retarded, causal" Green 
function that describes propagation away from (and not towards) a localized perturbation 
pulse. The restriction imposes a relation between the inlet conditions on A and d^A on the 
one hand and the bulk forcing / on the other hand. Thus, for example the inlet derivative 
dxA is fixed in the physically relevant solution of the second order (in space) CGLE by A 
at X = and the bulk forcing /. 

We have investigated the effect of arbitrary statistically stationary inlet fluctuations that 
are independent from the bulk stochastic forces and that are, e. g., advected through the 
inlet at X = into the system x > with the through-flow from an upstream subsystem, 
X < 0. Then we have considered a simple model for such an upstream subsystem that 
contains thermal fluctuations only: It is held at some subcritical driving /i < and the 
thermally generated fluctuations that are advected towards the inlet are described by the 
stochastic CGLE. Thus, fi allows to tune the strength of the advected inlet fluctuation 
spectrum. 
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If the driving fj, in the downstream part a; > is subcritical then all frequency modes 
are damped. Then the correlation spectrum C{x,uj) of amplitude fluctuations is at large 
downstream distances from the inlet no longer influenced by the inlet fluctuations at a; = 
and C{x, uj) approaches the translational invariant spectrum of an infinite system. However, 
for supercritical downstream driving n > inlet-forced as well as bulk-forced amplitude 
fluctuations arc growing in downstream direction for frequencies u within a finite band. 
Typically the contribution Cb{x,uj) caused by the additive bulk forces to the correlation 
spectrum C {x, lu) of amplitude fluctuations at the downstream location x rapidly outgrows 
with increasing x the contribution Cin{x,u!) that follows from the upstream model for inlet 
fluctuations. In the large-x regime far from the inlet both contributions, Cin{x,u!) as well 
as Cb{x,u), grow ~ ^-'^ im Ki{u))x ^^]-^ —JmKi > having a maximum close to u; = 0. But 
the prefactors in Cb and Cm are such that there the amplitude fluctuation spectrum C{x, u) 
typically is dominated by the bulk contribution Cb{x,u!). However, we have shown with our 
model how to realize experimental setups with subcritical upstream driving where the inlet 
fluctuations could be well controlled and well defined — uninfluenced by e.g. the pumping 
machinery — and tuned to be large enough to dominate over the bulk part Cb- 

Finally, a detailed comparison with a recent approach of SBH showed that the boundary 
condition that they use at the inlet implies fluctuations there which arc strongly correlated 
with the bulk forcing. This enhances the amplitude fluctuations of the growing modes and 
thereby tends to overestimate the effect of thermal bulk forcing. 
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APPENDIX A: GREEN FUNCTION 



We are interested in the Green function G{x, xq; t — to) of the GGLE as a solution of 

To{dt + vdx) - (1 + ico)Ai - (1 + ici)^ldl G{x, xq; t - to) ^ d{x - xo)S{t - to) (Al) 

subject to an arbitrary boundary condition G{x — 0, Xo] t — to) at the inlet. Because of the 
semiinfinite geometry G depends on x and the location xo of the (5-pulse separately. The 
temporal Fourier transform of G 

G{x, xo;u;) = H dte'''^^~*°'^G{x, xq; t - to) (A2) 

J —oo 

is conveniently written as the sum of two contributions 

G{x, Xo] uj) = e^^^'^G^x = 0, xo] a;) + G{x, xo] uj). (A3) 



Then the first term propagates the inlet condition in downstream direction and the boundary 
condition on the bulk contribution 
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:i+tc,)eo{Ki-K2) 



G(x,Xo;uj) 

e{x - xo)e'^'^'^-^"'> + e{xo - a;)e*^2(^-^«) - e^^-^^^-^^xo)! (^4) 



is G{x = 0,xo;u!) = 0. The physical solution ( p.l3|) of the linear stochastic CGLE can then 
be written in the form 



A(x,uj) = / dxoG(x,xo;uj)f(xo,uj) (A5) 

JO 

and 

POO pt 

A{x,t) = / dxo / dtQG{x,Xo;t — tQ)f{xo,tQ). (A6) 



The Xo-integration covers the whole spatial domain xq > 0. Causality implies contributions 
to come only from times —00 < to <t prior to the observation time t. The inlet conditions 
on G is defined by the requirement that for any given f{xo,uj) the integral 

dXoG{x = 0, Xq] U!)f{Xo, uj) = A{x = 0, uj) (A7) 

yields the imposed inlet condition A{x = 0, u) that can be chosen arbitrarily and in particular 
independent of /. 

For the sake of comparison with the Green function Gsbh used by SBH (cf. Appendix 
1]) we mention that G ([A3|) can be rewritten as the following decomposition 



G{x, xo; uj) = e*-^^"^ [G{x = 0, xo; uj) - Gsbh{-xo] uj)] + Gsbh{.x - xq; uj). (A8) 

Here 

Gsbh{x - Xo; uj) = . ' — \eix - xo)e'^'^^-^°^ + e{xo - x)e^^^(^-^'')l (A9) 

is the Green function used by SBH and 

Gsbh{-Xo; uj) = J e-^K,^o (Aio) 

(1 +zCi)4g(Ki - K2) 

is its value at the inlet. Note that the " space- retarded, causal" nature of the Green function 
can readily be read off from (|A9|) . It describes spreading off — not contraction towards — the 
perturbation pulse at Xq: the first (second) exponential in ( |A9| ) accounts for the spreading in 
downstream (upstream) direction with complex wave number Ki {K2). Downstream growth 
occurs for frequencies with ImKi{uj) < while the spreading in upstream direction is always 
damped given that ImK2 < 0. 

Since the linear operator appearing in the CGLE (|3.2|), i.e., the square bracket in ( [AI| ) is 
oc {dx — iKi){dx — iK2) it is clear from (|A8|) that both, G and Gsbh, solve Eq. (|Al|) , however 
with different inlet conditions. 
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APPENDIX B: COMPARISON WITH SBH 



1. Green function ^^^//(x — xo; cj) 
Here we derive the expression ( [A9D for 



GsBH{y; i^) 



with y = X — Xo by temporally Fourier transforming the Green function 

dk 



GsBH{y;r) 



^{l+iCo)/lT/TO 



To 



-oo 27r 



exp 



-(1 + ici)Qk^— + ik{y - vt) 
To 



(Bl) 



(B2) 



given in eqs. (C18, C19) of SBH for r = t — to > 0. Inverting the sequence of integrations 
one obtains 



GsBH{y]Uj) 



1 r°° dk 



To 
1 



°° dk 



.iky 



.iky 



1 



roo 

/ dr exp 
Jo 

POO 

/ dr exp 

^0 



T 



iujT + (1 + ico)^- (1 + ici)^Qk^ 

To To 



T 



ikvT 



[l + ic,)eo{k-K^){k-K2) 



To 



dk 



oiky 



1 + zci)^2 27r (A; - K^){k - K^)' 



(B3) 



The above integral along the real A;-axis can be calculated as a contour integral in the 
complex fc-plane with the method of residues. For y > (?/ < 0) one has to close the 
integration contour in the upper (lower) complex halfplane. 

The pole at K2 contributes only for y < since it lies with ImK2 < always in the 
lower half of the complex fc-plane. It yields the second term in ( pi|) . The pole at Ki lies 
for subcritical yU < in the upper half of the complex fc-plane. It causes the term oc e*^^^ 
for y > in ( pi| ) coming from the contour closed in the upper half of the complex fc-plane. 
However, when /i > this pole crosses the real fc-axis and moves into the lower /c-plane for 
frequencies within the band (|3.5| ) of growing modes. In order to retain its contribution in 
( pTD for 2/ > one has to deform the original integration path along the real fc-axis into a 
contour reaching into the lower fc-plane that encloses Ki such that this pole remains in the 
interior. 

Alternatively one can first perform the k-integration in ( P^ ) and then Fourier transform 
the result 



GsBH{y]^^) 



1 



2Co\/(l + «Ci)7rro 



°° dr 



exp 



iujT + (1 + ^Co)^i 



T_ _ {y - vrfro 

To A{l + lCi)eoT 



(B4) 



Using a/t as integration variable one finds with Eq. (3.472-1) of Gradshteyn and Ryzhik 
that 



GsBH{y\^) 



1 

— ( 

7 



(B5a) 



with 
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a 



VTo 



\ 



iuoTQ + (1 + ico)// 



(B5b) 

(B5c) 
(B5d) 



7 = 2/5(l + zci)eo^ 

Expressing a, /?, and 7 in terms of the eigenvalues K12 ( p^) and using the fact that \p\p jy 



sign{y) one verifies that (p5| ) agrees with (pi 



2. Amplitude fluctuations j45b//(x, w) 

The amphtude fluctuations of SBH are expressed with the help of the Green function 

GsBH as 



A 



SBH 



[x,t) 



00 i-t 

dxQ / dtoG 
J-00 



O^SBH 



{x -Xo;t- to)f{xo,to). 



(B6) 



Here we deal with the statistically stationary situation where the stochastic forces have been 
operating since the time to = —00. On the other hand, SBH enforce the initial condition 
AsBH{x,t = 0) = at time t = by integrating in ( |B6|) only from to = to to = t — cf. 
Eq. (C.20) in [^. Since we are interested however in comparing the longtime properties of 
their results with ours this difference does not play a role — the formulas labelled with the 
subscript SBH in this and the next subsection refer to the longtime limit of the SBH results 
in 0. 



In frequency space the longtime amplitude fluctuations of SBH follow with ( pi|) to be 



Asbh{x,uj) 



+ 



(B7) 



Note that the SBH solution is a very special one in the sense that it assumes the longtime 
amplitude fluctuations at the inlet to be of the form 



{l + zc,)eo{K,-K2) Jo 
{l + ^c^)eo{K,-K2y 



dxoe-''''-"fixo,u;) 



(B8) 



Thus, the inlet boundary condition is not treated as an independent quantity that can be 
imposed on the solution of the CGLE. Rather the inlet boundary value for A is related to 
and determined explicitly by the bulk forcing f{xo,uj). However, when we impose in our 
general solution (|3.14| ) the special inlet condition Asbh{x = 0,u) of SBH do we get the 



same result as SBH. In fact, our general solution (|3.14|) can be written in the form 

A{x, u) = AsBHix, iu) + [A{x = 0,uj)- AsBnix = 0, u)] (B9) 
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with an inlet boundary condition A{x = 0, u) that is still free to be chosen. 

Finally we mention that the derivative of the special SBH solution at the inlet 

d^AsBHix,uj) 1^=0 = , . ..2,^. -TT / dxoe~'^^'"'f{xo,uj) 

(1 + ici)^o{Ki - K2) Jo 

= iK2AsBHix = 0,u) (BIO) 
is such that also the SBH solution obeys the relation (|3.12|) 

i{d,-zK,)AsBH{x,u;) 1,^0+77^^^ = 0. (Bll) 

discussed in Sec. [Ill A| which prevents unphysical behavior for x ^ 00. 



3. Correlation function C5B//(x, cj) 



Using < Asbh{x,uj) [AsBnix^uj')]' 
spectrum one obtains from 



>= 27r6{uj—Lj')CsBH{x, u) to evaluate the correlation 



Csbh{x,uj) 



2(1 + cj)a\K,-K2 



ImKi 



(1- 



-2ImKix 



ImKo 



(B12) 



with inlet correlations being given by 



Csbh{x = 0,u;) 



-1 



K1-K2 P Jm/sTa 



(B13) 



These correlation spectra have to be compared with our result ( p.22|) and ( p.23|) . Note that 



F = 2^qTqF^^^ as mentioned already in Sec. |11 C| . 

To that end we show in Fig.^ amplitude fluctuation spectra at the inlet, x = 0, and 
further downstream, x = lO^o, for Re=2 and /i = ficonvf^ — 0.03. Full lines show our 
results using inlet fluctuation spectra (Sec. [IV A| ) obtained for different subcritical upstream 
driving fi as indicated. The most important difference is that in SBH the inlet fluctuations 
are determined by the bulk forcing / with a spectrum ( pi3[) that is fixed by F, Re, and 
H whereas in our more general case the inlet boundary condition on the amplitude A is 
still open. In the case that upstream fluctuations outside the main system generate our 
A{x = 0,uj) the latter is statistically independent of amplitude fluctuations Ab{x,uj) ( p.l4c[) 
that are generated in the bulk by the random forces /, i.e., < A{x = 0,u!) [Afy{x,uj')]* >= 

(PI)- 

This different correlation behavior influences not only the amplitude fluctuations close 
to the inlet but also the large-x correlation spectrum of the growing frequency modes with 
ImKi < 0. Only the correlation spectra of the damped modes with ImKi > are in both 
approaches the same for large x: Csbh{x 00, uj) = C{x oo,uj) = Coo(cu). ([3.26[) holds 
for the damped frequencies with ImKi > since at a; — > cxd these modes do no longer 
feel any influence of the inlet boundary condition. However, for the growing modes with 
ImKi < 
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Csbh{x 



2(i+c?)a 



-2ImKxx 



K1-K2 \^ImKi 



(B14) 



differs from 



C(x — 00, oj) 



D(uj) 



2(1 + cf)^^ I Kt - K2 p \ImKi ImK2 



-2ImK\x 



(B15) 



Thus, if one just replaces in ( [B15|) our inlet spectrum D{ijj) by Csbh{.x = 0,uj) ( [B13|) 
one does not recover the SBH result (|B14| ) since the SBH spectrum was derived for inlet 
fluctuations Asbh{x = 0,0;) ( pSD which are correlated with the bulk forcing / while our 
inlet boundary condition on A is chosen to be externally determined and thus independent 
of the bulk forcing. 

These additional correlations between inlet, x = 0, and downstream positions, x > 0, 
that are present in the SBH solution enhance the SBH spectra compared to our spectra. 
Consider, e.g., the case /i = —0.05 in Fig.^ for which our inlet fluctuation spectrum is 
even larger than Csbh{x = 0,u!). Nevertheless, the downstream amplitude fluctuation 
spectrum Csbh{x,uj) rapidly outgrows our C{x,uj). The effect of this correlation enhanced 
downstream growth of SBH can also be seen in Fig.^ — compare the dashed line for the 
peak height Csbh{x,u = 0) with our result C{x,uj = 0) for fi = —0.05. 

This effect can be measured quantitatively, e.g., by comparing the prefactors of the 
large-x exponential growth behaviour of (P14|) with that of (|B15| ). To that end we ignore 
the imaginary coefficients cq and ci in view of their smallness ||3^. Then the ratio of the 
spectral peak heights at = has the form 



where 



Csbh{x 00, lj = 0) 



C{x 



00, u 



0) 



2(1-^) ji-fi-fi/fi 



(BI6) 



< 



(B17) 



denote reduced control parameters. For not too large fi the ratio (pi6|) is close to 1. However, 
when fi approaches 1 it becomes large. The coefficients cq, ci remove the singularity present 
at /i = 1 in ( |B16| ). In Fig. 5 the ratio Csbh{x,lu = 0)/C{x,uj = 0) of the peak heigths of 
the full correlation spectra is shown at a; = 20^o versus reduced downstream driving ft. This 
ratio differs only slightly from the large- a; asymptote ( |B16| ) because of a small contribution 
from the last term in ( |3.22c|) . 

The inlet fluctuations resulting from the upstream part of the system would have to 
be very large, i.e., /x would have to be very close to zero in order to cause downstream 
fluctuation amplitudes that are as large as those of the SBH approach. Furthermore, if 
thermal equilibrium fluctuations were causing the inlet fluctuations as described by the 
unforced case /x = — 1 in our upstream model the discrepancy would be even larger. Thus, 
in agreement with an assessment made in the introduction of |2^ we conclude that the SBH 
result overestimates the effect of thermal noise mainly because the inlet boundary is not 
treated adequately. 
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FIGURES 

31) of amplitude fluctuations at the inlet normalized by TQD{t = 0) 



FIG. 1. Spectrum D{uj) 
versus reduced frequency lotq for several /i. Parameters [39| are for Re=2. 



FIG. 2. Downstream evolution of the reduced contributions Cin{x,uj ) ( ^32b| ) and Cb(x,u;) 
( 3.22c| ) to amplitude fluctuation spectra. Being practically symmetric in u> for our parameters 
Cin(,o/F (left side) and Ch^o/F (right side) can be continued symmetrically to the other side. Inlet 
fluctuations Cin{x = 0,a;) = D{ijj) are taken from our model with subcritical upstream driving 
/i = —0.1. Parameters are Re=2, ^ = = 0.03. 



FIG. 3. Downstream growth of spectral peaks. They are located very close to a; = for our 
parameters Re=2, /i = [f^^^^jl = 0.03. Thin lines show Cin{x,uj = 0)S,o/F coming from inlet 
fluctuations (Sec. |IV A[ ) obtained for different subcritical upstream driving For the SBH result 
Csbh{x^u) = Q)io/F (|B12 ) represented by the thick dashed line the inlet fluctuations are fixed by 
construction via /i, Re (cf. Appendix I 



FIG. 4. Reduced correlation spectra of amplitude fluctuations at (a) inlet, x = 0, and (b) 
further downstream, x = lO^o- Parameters ||3^ are Re=2, // = {fconvl'^ — 0.03. Full lines show 
our results using inlet fluctuation spectra (Sec. IV A) obtained for different subcritical upstream 
driving /i. Inlet fluctuations entering into the SBH result (dashed lines) are fixed by construction 
via fi, Re (cf. Appendix [B 3]). 



FIG. 5. Ratio of spectral peak heights of amplitude fluctuations according to SBH and to our 
work. The correlation spectra have almost reached their large- a; asymptotic form at a; = 20.^0 • 
Our spectra were evaluated using inlet fluctuations (Sec. IV A| ) obtained for different subcritical 
upstream driving /i. Inlet fluctuations entering into the SBH result are fixed by construction via 
/i. Re (cf. Appendix | 



3|). Parameters [3^ are Re=2, 



/2 = 0.03. 
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